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We study the consequences of the gravitino decay into dark matter. We suppose that the lightest 
neutralino is the main component of dark matter. In our framework gravitino is heavy enough to 
decay before Big Bang Nucleosynthesis starts. We consider a model coming from a five dimensional 
supergravity compactified on S 1 /Z2 with gravity in the bulk and matter localized on tensionless 
branes at the orbifold fixed points. We require that the dark matter, which is produced thermally and 
in the decay of Kaluza-Klein modes of gravitino, has an abundance compatible with observation. We 
deduce from our model that there are curves of constraints between the size of the extra-dimension 
and the reheating temperature of the universe after inflation. 



INTRODUCTION 



o 
o 

The five dimensional picture of the Universe has attracted much interest in the framework of what is called brane 
world cosmology. Several models have been proposed. It has been shown [l[ that a stabilization of the radion field 
K-j ! provides the classical Friedman equation. Many papers (see for instance: [2]) have considered a supersymmetric version 
of brane world which embed brane world in the framework of superstrings theories. In the present work we choose 
to work in a set-up where we assume that the radion is stabilized in a five dimensional supergravity compactified on 
S 1 jZi where matter and gauge fields live on the branes and gravity in the bulk (0|, Q). Low energy supersymmetry 
provides a natural candidate for dark matter if R-Parity Q is conserved and solves the hierarchy problem. We have 
to point out that in the present model the extra-dimension does not play a role to solve the hierarchy problem like 
m ADD model (@, Q) or Randall Sundrum @|. It is supersymmetry which plays this role. So the size of this 
extra dimension is not constrained like in 0| and actually we are dealing with smaller extra-dimension than in those 
scenarios. 

As a general framework we can choose a scenario in which susy breaking is mediated to the observable sector partly 
f*"*. ■ by gravity anomaly and partly by Scherk-Schwarz mechanism [lOj . It is a scenario which can avoid the appearance of 
tachionic masses which are present in pure anomaly mediated scenario and problems due to pure gravity mediation 
scenario. This mechanism provides high masses for gravitino (it means that the mass of the gravitino can be above 
Oh! lOTeV). 

The gravitino field has Kaluza-Klein excitations modes due to the presence of the extra-dimension. We suppose that 
all gravitino modes are produced after inflation during the reheating period by scattering effects in the primordial 
thermal bath. We suppose that the lightest mode - i.e. the zero mode - is heavy enough to mainly decay before the 
Big Bang Nucleosynthesis starts and we calculate its corresponding mass. Actually this is naturally the case in some 
scenarios of susy breaking (anomaly mediation or mix between anomaly and Scherk-Schwarz mechanism). Gravitinos 
modes decay into supersymmetric particles and standard model ones. If R-parity is conserved, all gravitino modes 
will give at the end of their decay cascade a Lightest Supersymmetric Particle (LSP) which is assumed to be the 
lightest neutralino. But not all of these decay products of gravitino will contribute to the relic density of neutralinos 
which is assumed to be the dark matter. Actually only the gravitino modes which decay after the thermal decoupling 
of neutralino contribute to the dark matter density. If a gravitino mode decays before the thermal decoupling of 
neutralinos it does not increase the number density of neutralinos since these gravitinos produce neutralinos which 
are in thermal equilibrium. So a finite number of gravitinos modes contribute to the non thermally produced dark 
matter. The total of the thermally produced and non thermally produced amount of dark matter is constrained by the 
evaluation of the dark matter content of the universe. As a consequence, we can draw curves of constraints between 
the size of the extra-dimension and the reheating temperature because the number of KK gravitinos modes is related 
to the size of the extra-dimension and the number density of gravitinos is related to the reheating temperature. We 
have chosen a reheating temperature in the range 10 5 GeV to 10 10 GeV. This range is quite natural for scenarios 
which allow baryogenesis through leptogenesis. There is another constraint on the size of the extra-dimension coming 
from the fact that there are also KK gravitons which can disturb BBN if the number of KK modes is too high. 
In the present work we first present the interactions between KK gravitino and MSSM, and discuss abundance and 
lifetime of gravitino. After this, we estimate the thermal production of neutralinos function of Xf = ^yr 2 - where Tf 
is the neutralinos freezing temperature and m/ sp the neutralino mass. We then present the model, the equations of 
constraints and check that KK gravitons are not a problem for BBN for Rr 1 > 1 TeV. Finally we show and comment 
the results. 
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II. INTERACTIONS BETWEEN GRAVITINO AND MSSM 

The four dimensional (reduced) Planck mass M is related to the five dimensional one M5 by the relation M 2 = 
ttRMq where R = jg- is the physical radius, r is an undetermined real constant and M 5 = . 
As matter and gauge fields live on the brane, those fields do not have a dependence on the fifth coordinate x 5 . The 
two gravitino fields (odd and even under Z%) have the following Fourier expansion: 
for the even fields under Z 2 : 



^™(x x ) + V2 £ C^V) cos(nM 5 x 5 ) 



(1) 



and the odd fields: 



We can define: 



r, dd (x\x 5 ) = ^^C d M d (^)sinMf 5a; 5 ), (2) 



7rr 

n=l 



^(x\x 5 ) = r; cn (x\x 5 ) + ^ dd (x\x 5 ) (3) 

,( /'n™ n ( a;A ) an d V'n d / d ( a;A ) are the Kaluza-Klein modes. 

Here are the relevant terms of interactions between matter, gauge fields and gravitino in the four dimensional 
supergravity lagrangian: 

~iu e ^^ VP ^\a) +^"^A (S) ) FW, (4) 

Where 4> are scalar fields, \ are chiral fermions, A are gauge fermions (gauginos) , F$ is the field strength tensor 

for the gauge boson A^. Indices i, j.. represent species of chiral multiplets and (a), (b).... are indices for adjoint 

representation of gauge group, e is the vierbein. gij* is the Kahler metric. 

This lagrangian Q describes the interaction of the mode with matter and gauge fields. 

We obtain the four dimensional lagrangian by integrating over the fifth dimension the five dimensional one. As the 
gauge fields and matter fields of the observable sector live on one brane (for instance the one located at x 5 = 0) , the 
part of the action which contains interaction of the gravitino field with the matter and gauge fields is : 

S = K Jd 4 xJ dx 6 S(x s )(-^=eg ij »V v ^x i <^o u M x>k " ^eg^.V^x 3 ^{x x , x 5 ) 

- l -e (Vv(x A , x 5 K^A (a) + ^(x A , x 5 )a v »a»\ (a) ) F<#) (5) 

Using the definition and inserting the Fourier expansion of gravitino {J) and Q into the above action, we 
find the interactions of each KK mode with matter and gauge fields after integrating on the fifth dimension. After 
redefining KK modes of gravitino ([3j, [4j): 

/even i ./.odd 

and 

<Po,, = i/$T (7) 



we find: 
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n=0 



- Wl e MV^,, + ^>^A (a) ) ^) (8) 

So we get an infinite sum of the Lagrangian Q . Each Kaluza-Klein mode has the same interaction with matter 
and gauge fields. This is not a surprise in the sense that graviton and gravitino are in the same supermultiplet and 
each KK gravitons has the same interaction with matter and gauge fields fill ]. 

Gravitino Kaluza-Klein modes differ only by their masses. They get a mass through the super-higgs mechanism. 
Other redefinitions of the gravitinos fields have to be performed to obtain the gravitino mass matrix. Those redefi- 
nitions do not play a role in the interaction part. After diagonalization of the mass matrix there are no interactions 
between the different KK modes. 

The masses of the modes n are related to the mass of the zero mode by the relation (01,3): 

M n = M + - (9) 
K 



III. ABUNDANCES 



In the present work we suppose[13| that gravitinos are produced during the reheating period after inflation by 
scattering processes [28] . All the MSSM particles are in thermal equilibrium except the gravitino which decouples as 
soon as it is produced because of its extremely weak interaction (gravitational one) . Ten reactions are considered to 
calculate the gravitino production. Infrared divergences appear in the calculation due to the exchange of massless 
gauge bosons in the t-channel. This problem is avoided by the method developed in [T3|. Another group ([l7].[18]) 
has recently done the calculation again with a slightly difference in the result. Abundance is given in [iji] taking into 
account production during the inflaton-dominated epoch and [18] has the same result without taking into account 
production during inflation. This result is given for masses of gravitino much higher than gauginos masses but the 
calculation of the creation term in the Boltzmann equation is made with particles without masses: their mass is 
supposed negligible compared to the average energy in the center of mass of each reaction. The formula is : 



Y, 



3/2 



1.9 x 10" 



Ik 



10 10 GeV 



1 + 0.045 In 



10 10 GcV 



1 - 0.028 In 



10 10 GeV 



(10) 

where Y 3 / 2 = —j^-, "3/2 ls the number density, s is the entropy density and Tr the reheating temperature. The 
quantity Y — j is the density per comoving volume. 



We have to take into account in our calculation the different masses of gravitino modes. We find as a good 
approximation this rule for the abundance of the different modes: 



Y 3/2 = Y 3/2 > for Mk < Ik and 

Y 3 ) 2 = , for M k > T R (11) 

where k represents the index of the KK mode, M k is the mass of the k th gravitino mode and Y 3 , 2 its abundance. 
In this model all particles are considered massless compared to the reheating temperature scale except the gravitino. 
All particles in the thermal bath are relativistic at Tr except the gravitino modes which can be non relativistic if 
heavy enough. 



IV. DECAY OF THE GRAVITINO MODES 



In a previous section we showed that each gravitino mode had the same interaction with matter and gauge fields. 
The only difference arises from the mass of the different modes. 
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Gravitino interacts gravitationally with MSSM particles: its life time is long compared to the one of other particles. 
We have calculated the life time for heavy gravitino (masses > 10 TeV) using [TjJ: 



^ = ^10 7 x(^)"sec (12) 

In our model R-parity is conserved and the LSP is the lightest neutralino. 
Several decay channels are possible. Heavy gravitino decays mainly into squark-quark pairs or gluon-gluino pairs 
We checked that finally one gravitino produces one non-thermal neutralino and not less because number density 
of susy particles produced by gravitino decay is so weak (Y » 10~ 12 ) that they will more likely decay than interact 
to produce standard model particles: in the Boltzmann equation the decay term is proportional to Y and the 
annihilation term is proportional to Y 2 . We also checked that annihilation of the LSP can be regarded as negligible 
because it will concern only the LSP produced by one gravitino mode when the total of modes is three. Gravitino 
does not produce more than one LSP because even high energy quark or gluon produced by gravitino decay if they 
interact between them or with a quark or a gluon of the plasma will almost all time produce quarks or gluons rather 
than neutralinos as processes with strong interactions are favored compared to the ones with weak interactions. We 
also checked that the decay of gravitino modes does not increase entropy density because the energy density of each 
gravitino mode at the time of its decay is never higher than the radiation energy density at the same time: it is 
because all the gravitino modes decay before T — 1 MeV. 



V. NEUTRALINOS 

In our model the LSP is the lightest neutralino. We choose to work with a mass of LSP equal to 120 GeV. This 
analysis can be easily rescaled for another choice for the mass: we also show results for a mass of the LSP equal to 
200 GeV in the appendix. The dark matter density is [23] : 

0.106 < Q h 2 < 0.123, (13) 
with a central value around 0.114 

In the graph representing evolution of Yi sp as a function of the quantity x — ^J 2 - there are two zones. The first 
zone is the equilibrium zone in which neutralinos are in thermal equilibrium and in which Y decreases. The second 
zone is the zone after the freeze-out of neutralinos and in which Y is constant. The freeze-out occurs at Xf. We call 
flth the thermal density of neutralinos. We find this approximate relation between flth and xj : 

Q th h 2 = 3.61 10 6 ^ x 2 e~*t (14) 

To establish this relation we used [21]. This corresponds to the relations below in which we have neglected the 
double log part. 

x f = ln[0.038(rc + l)-j^m p z m ;sp a ] - (n + -) ln[ln[0.038(n + l)-^-m p/ m lsp a ]] (15) 
9* 1 9* 



and 



(n + l)x n+1 

n th h 2 = 1.07 x 10 9 UT ^ , (16) 

(g*s/g* )m p t cr 

with n = 1 (for p — wave annihilation), g* = g*g = 90, g = 2 

where g± is the number of effectively massless degrees of freedom, m p i is the Planck mass. 

We choose different values of fl t h h 2 and complete with the non thermal production coming from the gravitino 
decay. We call this non thermal production Aft h 2 . 

0.106 < n t h h 2 + AQ h 2 < 0.123 (17) 
As one gravitino produces one neutralino, we can write: 
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FIG. 1: flth h 2 function of Xf. Q,dm h 2 is between the two straight lines. 
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An h 2 = 



Pc 



Yl F 3/2 



(18) 



k=0 



The index n corresponds to the last mode to be taken into account. It is the mode decaying just when the LSP 
decouples from the thermal bath. Only the gravitino modes decaying after the thermal decoupling of the neutralinos 
contribute to increase the quantity of neutralinos. 



The thermal density Q t h h 2 and x / are related by the equation lfl4|) . Knowing the thermal density we can calculate 
the non thermal density to respect the requirement (fTF)) . 

Knowing i/, we can calculate Tf and then the mass of the last gravitino mode taken into account. 

To do so we need a relation which allows us to convert temperature into time. In a radiation dominated Universe 

assuming that the entropy per comoving volume is conserved the relation between t and T is: 



Where M is the reduced Planck mass. This result (fl9| is obtained by solving the Friedman equation in a radiation 
dominated Universe assuming constant entropy per comoving volume. 

We use relation lfl9|) to calculate the mass of the mode n but also to calculate the mass of the mode. For the mode n, 
(J* = 90 and for the mode, g* = 10. We write equality between the life time of the gravitino mode lfT2|) and the age 
of the universe[29] fl9|) (after having converted (fT9| into second) and we can calculate the mass of the corresponding 
gravitino. 

The gravitino which mainly decays at Tf corresponds to the last mode to take into account in our Kaluza-Klein tower. 
We call it the mode n. 

The mass of the zero mode is calculated the same way : we just put as a condition that it mainly decays at T = 1 
MeV not to disturb BBN. For a temperature equal to 1 MeV we find with (fl9|) t = 0.76 s. We then calculate M with 
(p| by writing: r = 0.76 s. We get M = 26.410 TeV. 

The mass of the mode n depends on the choice of Xf which depends on the thermal density flth h 2 . 



VI. MODEL 



A. 



Masses 




(19) 



B. Numerical results 

We chose three different values of the density Q t h h 2 . 
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The first chosen value corresponds to the central value h 2 = 0.114. This value is the value corresponding to 
our maximal thermal contribution. 

The second value is calculated by choosing the maximum value for cr . This cross section is generically bounded [HJ 
by: 

vo < (20) 

m Lsp 2 

With a ~ 10~ 2 . By using the maximum value for the annihilation cross section we find our minimal thermal 
contribution. 

The last value is a medium value: we chose an Xf value between our two Xf extrema and calculate the thermal 
density corresponding to this x/ value. 
Below is the table of values: 



Cases 


x f 


T>(GeV) 


n th h 2 






M ft (GeV) 


Case 1 


28.78 


4.17 


0.114 




0.009 


9.84 x 10 6 


Case 2 


29.60 


4.05 


0.053 


0.053 


0.070 


9.66 x 10 6 


Case 3 


30.42 


3.94 


0.025 


0.081 


0.098 


9.48 x 10 6 



TABLE I: The three numerical cases for mi ap = 120 GeV 



1. Equations of constraints 

Using equation (fl8|) and replacing Yy 2 by its value ifTTj) . we obtain the non-thermal production of LSPs. 
We have to distinguish two cases. 

The first case is when the mass of the last mode to be taken into account (so the mode that we have called n and for 
which we have calculated the mass in the last section) is less than Tr or equal to Tr so Tr > M n . In that case if we 

! 3/2 ~~ x 3/2 



look at the equations (flT|) we deduce that Y^, 2 = Y®, 2 and so we obtain from equation (fT8|) : 



mi sp (n + l)Y° /2 s h 2 

An h 2 = — — (21) 

Pc 



We then should replace n by its value function of R , M™ and M°: 

M fl - M° 



(22) 



Let us keep in mind that n is an integer number which means that R 1 which is one of our variable must have a value 
which allows the ratio R ~i to be an integer. 

So we will replace n by its value after having isolating it in one side of equation (f2~Ij) . We define a function I which 
acts on real number and extracts the integer part. 
We obtain: 



R- 1 = (Af fi - M°) x 



Aft h 2 



misp s h 2 1.9 x 10- 12 x ^ [1 + 0.045 In (^)] [l - 0.028 In (^)] 



(23) 



In the above equation (|23)) we use Table [J to fix the limit values on AO h 2 . 
The above equation is valid for Tr > M n and for a number of modes at least equal to 2. This last condition implies 
a limit on Tr: above the temperature corresponding to a number of modes equal to 2 there is only one mode which 
produces all the non thermal neutralinos. We calculated those two limits on Tr for the three cases. The results are 
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case 1 


case 2 


case 3 


1.45 10 9 GeV 


1.09 10 10 GeV 


1.52 10 10 GeV 



TABLE II: Maximum reheating temperature 



presented with the graphs. 

We also calculated the maximum reheating temperature allowed if there is only one mode for the three cases. The 
results are shown in Table [ill and converge to the limit found in [53] in the limit of vanishing thermal relic density. 

If T R < M n , it means that the mass of the last mode that should be taken into account has an abundance equal 
to zero using the rule given in (TTJ) . In that case the last mass that it is effectively taken into account is equal to Tr. 
We obtain the result for this case by just replacing M n by Tr in equation f23|) . 



VII. GRAVITONS 



As KK gravitinos are produced, there are also KK gravitons. We have to check that they do not disturb BBN 
since they are heavier than in the case of large extra dimensions usually considered. In our framework the usual 
astrophysical and cosmological constraints disappear like in the case of [24j]. But there is a new problem coming from 
the fact that BBN should not be disturbed. The mass of the KK graviton is given by: 



™fc = ^ (24) 

So the mode 1 is the first massive mode with a mass equal to j=r. In the following we check that for > 1 TeV, 
BBN is not disturbed. For larger R the presence of KK gravitons starts affecting BBN as shown with our approximate 
method. A more accurate study is needed to check what is precisely the maximum size which can be allowed for the 
radius. This calculation is beyond the scope of the present work and will be discussed elsewhere. 
The production equation of KK gravitons is given by [9j and [25J]. We can use the same equation: 

llm 5 T 

S¥m = 128^p (25) 

Where K\ is the modified Bessel function of the first kind, m is the mass of the graviton, T the temperature and 
M the reduced Planck mass. 

By integrating this equation between the reheating temperature and temperature T below 1 MeV, we find: 

oo 

Y m (T) = , 1485 =-7= m [ x 3 K 1 (x)dx (26) 

v ' 256tt 5 M g*V 2 g* s J 

m/T R 

where g* and g^ are taken constant equal to 10 since most of the lifetime of the gravitons that we consider is after 
T = 1 MeV. 

We calculate the lifetime of KK gravitons with a mass above 1 TeV using the partial decay rates given in (nj. We 
find : 

t = 3.310 x GeV" 1 (27) 

Ttl 

We then calculate the mass of the graviton which decays at T = 1 MeV using equation lfl9|) and equation (|27")l . 
This graviton is the last one to be taken into account in our check. We find m — 37.4 TeV and we calculate its 
lifetime r = 0.75 s. As the radius is equal to 1 TeV the last mode that we consider in this study has a mass equal to 
38 TeV. 

To check if KK gravitons disturb or not BBN we use the curves given by Jedamzik in [26] . Those curves give a limit 
on the density of a massive decaying particle (if it did not decay) as a function of its lifetime. We have calculated the 
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density of gravitons at our epoch if they did not decay. We also have to know the hadronic branching ratio. We find 
Bh = 0.70 for gravitons with mass above 1 TeV using the partial decay rates given in fTTj |. 

Those curves are made for one particle decaying and not like in the present case for a tower of particles. We can 
however roughly divide the study in two zones : one before r = 100 s and the other after t = 100 s. 
In the first zone (i.e r > 100 s) we take as a limit this estimate value: fl h 2 = 5 10 -5 and before 100 s we take as a 
limit n h 2 = lCT 1 . 

The mass of the graviton whose lifetime is 100 s is 7.33 TeV. So we sum the density contribution of the first seven 
modes of the first zone and we check that this sum is less than 5 10~ 5 . We also sum the density contributions of the 
31 modes of the second zone and check that this sum is less than 10 _1 . 



The density for gravitons is : 



k—n 

E ™fc Y k sq h 2 k=n oo 

nh 2 = ^i = s AJLy m i 1485 / x s K{x)dx {28) 

Pc Pc fr{ fe 256 7r 5 M 5 *V2 5 * J 1 > 

m k /T R 

Where we have used equation (|26| in which we have replaced m by rat and Y m by Yk since each mass is associated 
to a mode. We know that: 



oo 

/ <«1 (29) 



m k /T R 

This integral is equal to 4.71 if m = or if m « Tr, If m — Tr, the integral is equal to 4.47. In the range of 
values for m and Tr that we have, we can choose to set the value of the integral to 4.71: for masses below 7 TeV, it 
is obvious that iti/Tr « 1 and for masses above 7 TeV and close to 38 TeV, the error that we make increases the 
density. 

This implies : 



so h 2 5pJ 1485 m 



A Vps 4.71— > = 7-77; — (30) 

p c 256 ir>Mg* l l 2 g* s y ' 

We can replace rri]~ by 4 and after doing the sum over k we find: 

h 2 « 5.52 10 -i4"d + "Kl + 2n) (31) 

Note that n = m n R and therefore fl h 2 is proportional to R. So the larger is the radius, the larger is fl h 2 . 
The mass of the last mode of the first zone is 7 TeV so after having replaced n by m n R we can evaluate the contribution 
of the first seven modes. We find: 

Q h 2 = 4.64 10" 5 < 5 10~ 5 (32) 

The result is below the limit even if it is rather close to it. But this limit is also not completely defined for our 
particular study in which there is a tower of massive particles. We conclude as a first estimate that the density of the 
seven first modes of gravitons does not disturb BBN. 



The mass of the last mode of the second zone is 38 TeV so after having replaced n by m n R we can evaluate the 
contribution of the 31 modes. In our calculation we take into account the seven first modes whose contribution is 
negligible. We find: 



n h 2 = 6.30 10~ 3 < 10" 1 (33) 

The result is well below the limit. We conclude as a first estimate that the density of the 31 modes of gravitons 
does not disturb BBN. 

We can conclude this section by saying that the curves drawn from gravitinos constraints on dark matter density 
are valid for R^ 1 > 1 TeV. This new constraint comes from KK gravitons. A more precise study has to be done to 
check if R^ 1 > 1 TeV is a stringent limit or if it is possible to have a larger radius. 
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VIII. RESULTS 

In this section we present our results for the three different cases. Let us remind that the first case corresponds to 
the case in which il t h h 2 is our maximum 0.114 and so the non-thermal part has to be less than 0.009. The second 
case is our intermediate case in which Q t h h 2 = 0.053 and so the non-thermal part has to be more than 0.053 and 
less than 0.070. The third case corresponds to the minimum value for Ctth h 2 so 0.025 and the non-thermal part has 
to be more than 0.081 and less than 0.098. The range of values that we have chosen for Tr is from 10 5 to 10 10 GeV. 
All these cases are treated with m/ sp = 120 GeV. We also treat three cases presented in the appendix for mi sp = 200 
GeV. 

The first general comment on the curves is that they are all increasing. The lower is the reheating temperature 
the higher is the number of modes needed (and so the smaller is R~ x ) because the gravitino density (and so the 
non-thermal part of neutralino density) is nearly proportional to the reheating temperature. 

We can check that the curves of the first case are below the curves of the second case which are also below the curves 
of the third case: the more dark matter has to be produced the more modes is needed and the smaller is. We 
can also check that all curves are converging to nearly the same RT 1 : it is actually the maximum mass difference 
between the first mode which mainly decays at T = 1 MeV and the last mode which mainly decays at T/ and we have 
seen that our three different M n have close values (see Table HJ . We notice in all the cases that a discrete structure is 
visible above a different temperature for the three cases. Actually the discrete structure is always present if we zoom 
enough on the graph. For high scale temperatures, this structure is always visible on our graphs since the number of 
modes is low. 

We also trace on the figures the limit > 1 TeV to avoid perturbation of BBN by KK gravitons. 



R 1 (GeV) 



10000 : 



1000 , 



100 : 



10 r 




' . . ' . — . ■ ■ i 1 , . 1 . — . . Tr (GeV) 

100000. 200000. 500000. lxio 6 2.X 10 6 5.X 10 6 l.xlO 7 



FIG. 2: Case 1. Tr less than 9.84 10 6 GeV. The excluded zone is below the diagonal curve and below the straight line R 1 = 1 
TeV if the KK gravitons constraint is taken into account. 

The figure ^ is the case 1 for reheating temperatures below 9.84 10 6 GeV. This last temperature corresponds to 
the limit case M n = Tr. The excluded zone is below the curve. 

The maximum size for R corresponds to the minimum value for R , Without taking into account KK gravitons 
constraint, the minimum would be R^ 1 = 3.35 GeV for a reheating temperature equal to 10 5 GeV. So R should be 
equal or smaller than 5.89 10 -15 cm for Tr — 10 5 GeV. Taking into account the limit R^ 1 > 1 TeV, R has to be 
equal or smaller than 1.97 10~ 17 cm. For a reheating temperature of 9.84 10 6 GeV, R has to be equal or smaller than 
3.46 10~ 19 cm which corresponds to 5.71 10 4 GeV for H -1 . 

At the scale of the graph we cannot distinguish the discrete structure of the function. The maximum number of 
modes allowed should be equal to 21985 for Tr = 10 5 GeV without the constraint coming from KK gravitons and to 
172 for T R = 9.84 10 6 GeV. 
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FIG. 3: Case 1. 9.84 10 6 GeV < Tr < 4 10 7 GeV. The excluded zone is below the curve. 

The figure J3]) is the case 1 for reheating temperatures between 9.84 10 6 GeV and 4 10 7 GeV. The excluded zone is 
below the curve. For a reheating temperature of 4 10 7 GeV, R has to be equal or smaller than 7.84 10~ 20 cm (which 
corresponds to 2.52 10 5 GeV for We distinguish a discrete structure on the graph. The maximum number of 

modes allowed for Tr = 4 10 7 GeV is 39. 



IT 1 (GeV) 
1 . x 10 7 



Tr (GeV) 



FIG. 4: Case 1. 4 10 7 GeV <T R < 7.396 10 s GeV. The excluded zone is below the curve 



The figure (@]) is the case 1 for reheating temperatures between 4 10 7 GeV and 7.40 10 8 GeV. The excluded zone 
is below the curve. For a reheating temperature between 4.98 10 8 GeV and 7.40 10 s GeV, the maximum number of 
modes allowed is 2 and R has to be equal or smaller than 2.01 10~ 21 cm (which corresponds to 9.817 10 6 GeV for 
Above the temperature 7.40 10 8 GeV, only a single mode is allowed and R has to be smaller than 2.01 10~ 21 
cm. The maximum allowed reheating temperature is given in Table HU 1.45 10 9 GeV. Above this temperature the 
number of neutralinos is above the observational limit ([T4ll . 
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FIG. 5: Case 2. T R less than 9.66 10 6 GeV. Only the band between the two diagonal curves is allowed. The zone below the 
straight line -R -1 = 1 TeV is excluded if KK gravitons constraint is taken into account. 

The figure is the case 2 for reheating temperatures below 9.66 10 6 GeV. Only the band between the two 
diagonal curves is allowed. It means that R has to be large enough to provide the necessary amount of dark matter 
but not larger than the value corresponding to the maximum amount of dark matter. If KK gravitons constraint is 
taken into account, the zone below R^ 1 = 1 TeV is excluded. 

Without this constraint, should be between 0.43 GeV and 0.57 GeV for Tr = 10 5 GeV: it means that R should 
be between 3.47 10~ 14 cm and 4.58 10~ 14 cm. For T R = 9.66 10 6 GeV, R" 1 has to be between 7.02 10 3 GeV and 
9.28 10 3 GeV, so R has to be between 2.13 10~ 18 cm and 2.81 10~ 18 cm. For T R = 10 5 GeV, the number of modes 
should be between 129427 and 170957. For Tr — 9.66 10 6 GeV, the number of modes has to be between 1038 and 
1372. 

With KK gravitons constraint, the minimum value for is 1 TeV so the maximum value for R is 1.97 10~ cm. 
It implies that there is a minimum value for Tr corresponding to this value: the minimum value for Tr is 3.26 10 6 
GeV. Lower temperatures are excluded. 

The figure ((6|) is the case 2 for reheating temperatures between 9.66 10 6 GeV and 10 s GeV. Only the band between 
the two curves is allowed. For T R = 10 8 GeV, i?~ 4 has to be between 8.03 10 4 GeV and 1.06 10 5 GeV so R has to be 
between 1.86 10 -19 cm and 2.46 10~ 19 cm and the number of modes between 91 and 120 modes. 

The figure ([7|) is the case 2 for reheating temperatures between 10 s GeV and 4.19 10 9 GeV. Only the band between 
the two curves is allowed. For a reheating temperature between 3.70 10 9 GeV and 4.19 10 9 GeV, the number of 
modes has to be equal to 2 and R has to be equal to 2.05 10~ 21 cm (which corresponds to 9.634 10 6 GeV for R^ 1 ). 
Above the temperature 4.19 10 9 GeV , only a single mode is allowed and R has to be smaller than 2.05 10 -21 cm. 
The maximal allowed reheating temperature is given in Table Ull 1.09 10 10 GeV. Above this temperature the density 
of neutralinos is above the observational limit given in lfl4|) . 

The figure ((8]) is the case 3 for reheating temperatures below 9.48 10 6 GeV. Only the band between the two 
diagonal curves is allowed as in case 2. If KK gravitons constraint is taken into account, the zone below R^ 1 = 1 TeV 
is excluded. Without this constraint, should be between 0.31 GeV and 0.37 GeV for Tr — 10 5 GeV: it means 
that R should be between 5.33 10~ 14 cm and 6.44 10~ 14 cm. For T R = 9.48 10 6 GeV, i?~ 4 has to be between 4.81 10 3 
GeV and 5.82 10 3 GeV, so R has to be between 3.39 10" 18 and 4.10 lO" 18 cm. For T R = 10 5 GeV, the number of 
modes should be between 198804 and 240333. For Tr — 4.47 10 6 GeV, the number of modes has to be between 3568 
and 4314. 

With KK gravitons constraint, the minimum value for R~ x is 1 TeV so the maximum value for R is 1.97 10~ cm. 
It implies that there is a minimum value for Tr corresponding to this value: the minimum value for Tr is 4.02 10 6 
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FIG. 6: Case 2. 9.66 10 6 GeV <Tr< 10 s GeV. Only the band between the two curves is allowed. 
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1 

PIG. 7: Case 2. 10* GeV <Tr< 4.190 10 a GeV. Only the band between the two curves is allowed. 

GeV. Lower temperatures are excluded. 

The figures iJH) and (flOjl are the case 3 for reheating temperatures between 4.47 10 6 GeV and 10 8 GeV for the 
figure J9J and between 10 8 GeV and 2 10 s GeV for the figure l|10p . Only the band between the curves is allowed. For 
T R = 2 10 8 GeV, R- 1 has to be between 1.15 10 5 GeV and 1.39 10 15 GeV so R has to be between 1.42 1CT 19 cm and 
1.71 10 -19 cm and the number of modes between 68 and 82 modes. 

The figure ifTTj) is the case 3 for reheating temperatures between 2 10 8 GeV and 6.38 10 9 GeV. Only the band 
between the two curves is allowed. For a reheating temperature between 5.17 10 9 GeV and 6.38 10 9 GeV, the number 
of modes has to be equal to 2 and R has to be equal to 2.09 10~ 21 cm (which corresponds to 9.46 10 6 GeV for 
Above the temperature 6.38 10 9 GeV, only a single mode is allowed and R has to be smaller than 2.09 10~ 21 cm. The 
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FIG. 8: Case 3. Tr less than 9.48 10 GeV. Only the band between the two diagonal curves is allowed. The zone below the 
straight line = 1 TeV is excluded if KK gravitons constraint is taken into account. 
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FIG. 9: Case 3. 9.48 10 6 GeV <Tr< 10 s GeV. Only the band between the two curves is allowed. 



maximum allowed reheating temperature is given in Table ITU 1.52 10 10 GeV. Above this temperature the density of 
neutralinos is above the observational limit given in l[T4j) . 
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FIG. 10: Case 3. 10 s GeV <T R <2 10 s GeV. Only the band between the two curves is allowed. 




FIG. 11: Case 3. 2 10 8 GeV <T R < 6.385 10 9 GeV. Only the band between the two curves is allowed. 



IX. CONCLUSION 



We have shown in this paper that for cosmological models with high reheating temperature i.e. 10 5 GeV to 10 10 
GeV - in the framework of a 5D supergravity compactified on S 1 /Z2 where matter and gauge fields live on tensionless 
branes at the orbifold fixed points - there are curves of constraints between the size R of the extra-dimension and 
the reheating temperature. It comes from the assumption that dark matter is made by the lightest supersymmetric 
particle which is supposed to be the lightest neutralino and that neutralinos density is a sum of a thermal production 
and a non thermal production from gravitinos decay. Gravitinos in the model do not disturb BBN because they are 
heavy enough to decay before BBN starts. Heavy gravitinos are natural in a certain class of Susy breaking models. 
For instance models which avoid problems from pure anomaly mediation and pure gravity mediation. The framework 
of the model described in this paper can be linked to Horava-Witten M-theory where a 5 dimensional stage of the 
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universe appears in which bulk fields are gravitational and where supersymmetry is natural but also to theories of 
baryogenesis through leptogenesis which imply large reheating temperature. 

The results that we obtain are independent from the susy mass spectrum since the gravitino is heavy enough to make 
negligible the influence of other susy particles. 

The results show that the size of the radius R is not only bounded by a maximum value but also by a minimum 
value in a wide range of possible values for the thermal production of neutralinos and in a wide range of values for 
the reheating temperature. To obtain a minimum size for the radius is a new result. 

Kaluza-Klein modes of graviton are also present and may disturb BBN. We checked with an approximate method 
that for above 1 TeV it is not the case. This already implies a bound on the reheating temperature which can 
not be lower than a minimum value in the cases where the radius R is bounded by a minimum value. 
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APPENDIX: NUMERICAL RESULTS AND CURVES FOR m lap = 200 GEV 



Cases 


x l 


I>(GeV) 


fit,, h 2 






M„(GeV) 


Max.Reheat.Temp(GeV) 


Case 1 


29.32 


6.82 


0.114 




0.009 


1.37 x 10 7 


8.84 x 10 s 


Case 2 


29.6194 


6.7523 


0.087 


0.019 


0.036 


1.36 x 10 7 


3.45 x 10 9 


Case 3 


29.91 


6.69 


0.066 


0.040 


0.057 


1.35 x 10 6 


5.37 x 10 9 



TABLE III: The three numerical cases for 200 GeV 
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FIG. 12: Case 1 for m lsp = 200 GeV. T R less than 1.37 10 7 GeV. The excluded zone is below the curve. 

The figures (12]), (13|) and (JT4J) represent the case 1 but for m lsp = 200 GeV instead of 120 GeV. We notice that 
the curves of that case are above the curves for the case 1 with mi sp — 120 GeV. We can explain this by looking at 
equation (|23|) : we notice that is nearly proportional to mi sp . Physically it means that if the elementary energy 
of one LSP increases, less number density of LSPs is required and so can be bigger. M n , Tf and the maximum 
reheating temperature are slightly different from the case with mi sp — 120 GeV (see Table Hll|) . 

The figures of case 2 (figdUto fig fl8|) and case 3 (figHUto figETJ for mi sp = 200 GeV have the same general behavior 
as cases 2 and 3 for mi sp — 120 GeV: the curves are just above the ones for m/ sp = 120 GeV. 

For the case 2 the minimum reheating temperature is 1.55 10 6 GeV if the constraint on KK gravitons is taken into 
account (see fig fT5|) . For the case 3 the minimum reheating temperature is 2.216 10 6 GeV if the constraint on KK 
gravitons is taken into account (see fig fl9)l . 
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FIG. 15: Case 2 for mi sp — 200 GeV. Tr less than 1.36 10 7 GeV. Only the band between the two diagonal curves is allowed. 
The zone below the straight line R^ 1 = 1 TeV is excluded if KK gravitons constraint is taken into account. 
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FIG. 16: Case 2 for m !sp = 200 GeV. 1.36 10 7 GeV <T R < 2.5 10 7 GeV Only the band between the two curves is allowed. 
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FIG. 18: Case 2 for m Up = 200 GeV. 9.40 10 8 GeV <T R < 1.75 10 9 GeV The zone below the curve is excluded. 
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FIG. 19: Case 3 for mi sp = 200 GeV. Tr less than 1.35 10 7 GeV. Only the band between the two diagonal curves is allowed. 
The zone below the straight line R~ x = 1 TeV is excluded if KK gravitons constraint is taken into account. 
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FIG. 20: Case 3 for m !sp = 200 GeV. 1.35 10 7 GeV < T R < 5 10 7 GeV Only the band between the two curves is allowed. 



